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RECENT  RESEARCH  IN  COMPOSITE  AND  SANDWICH  PLATE  DYNAMICS 


C/W.  Bert* 


Abstract  - This  paper  surveys  literature  concerning 
dynamics  of  plate-type  structural  elements  of  either 
composite  materia I or  sandwich  construction.  Papers 
from  1976  through  early  1979  are  reviewed.  Special 
attention  is  given  to  rectangularly  orthotropic,  cylin- 
drically  orthotropic,  and  anisotropic  plates;  laminated 
plates;  thick  and  sandwich  plates;  and  nonlinaarities 
Free  vibration,  harmonic  and  random  forced  vibra- 
tion, thermally  and  flow  induced  vibration  (flutter), 
and  impact  are  also  treated. 

' \ 

The  fundamental*  of  the  mechanics  of  composite 
and  laminated  plates  have  been  discussed  in  a pre- 
vious survey  (1]  and  are  not  repeated.  No  books 
have  yet  been  published  in  the  subject  area:  the  sur- 
vey thus  consists Jp^ply  of  papers  in  the  open  litera- 
ture and  some  reports.  The  following  topics  are 
excluded:  strictly  in-plane  motion,  acoustic  wave 
transmission,  and  failure  due  to  impact  and  fatigue 
loadings. 


SMALL-DEFLECTION  MOTION  OF  THIN. 
SYMMETRICALLY  LAMINATED, 
RECTANGULARLY  ORTHOTROPIC  PANELS 


These  panels  are  the  simplest  and  most  extensively 
investigated  of  thin  composite-material  plates  and  are 
often  called  orthotropic  plates.  The  panels  can  be 
an  aligned  single  layer,  an  aligned  parallel-ply  lamin- 
ate, or  a symmetrical  cross-ply  laminate.  Laminate 
symmetry  about  the  plate  midplane  ensures  the 
absence  of  bending-stretching  coupling;  aligned 
orthotropic  material  symmetry  prevents  shear- 
normal  coupling.  Thus,  in  the  absence  of  applied 
shear  and  in-plane  loading,  the  governing  equation 
of  motion  can  be  written  as 


D»*w'xxxx+  2(Dia 

(1» 

+ Oj  , w Yyyy  + phw^  - p « 0 

The  Djj  are  the  plate  flexural  and  twisting  rigidities, 
h is  the  plate  thickness,  p is  the  normal  pressure, 
t is  time,  w is  the  plate  deflection,  ( ),xxyy  denotes 
34(  )/dx,3y,(  x and  y are  rectangular  position 
coordinates  coinciding  with  the  material-symmetry 
directions,  and  p is  the  mean  density. 

In  terms  of  energy  the  analogous  expression  is 

U,  + Ue-T  (2) 

The  internal  strain  energy  Uj,  work  Ue  done  by  the 
external  force,  and  kinetic  energy  T are  given 
respectively  as 

Uj  m (1/2)  //  (D»  jw  ,xx+  2D|^w,^yW^y 
+ Dj2w,,yy+  4D««wJ ,Xy)dx  dy 

(31 

Ue  * - //  p w dx  dy 
T -<1/2)//phw*ndx  dy 
Free  Vibration 

A method  for  deducing  orthotropic  solutions  from 
corresponding  isotropic  solutions  was  introduced 
independently  in  Japan  [2, 3]  and  in  the  U.S.  [1,4]. 
This  approach  enables  the  vast  literature  on  free 
vibration  of  isotropic  plates  of  various  planforms  to 
be  utilized  to  predict  approximate  natural  frequencies 
for  orthotropic  plates  of  the  same  planforms.  The 
accuracy  of  the  method  has  been  verified  for  the 
specific  geometries  listed  in  the  Table.  The  effects 
of  biaxial  inplane  foods  (tensile  or  compressive) 
were  later  incorporated  into  the  method  [5] , thus 
generalizing  previous  work  on  buckling  of  orthotropic 
plates  [6, 7] . 
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Table.  Applications  of  Isatropic-to-Orthotropic 
Deduction  Technique 


Planform 

Boundary 

Reference 

Geometry 

Conditions 

Number 

Rectangular 

Fifteen  Combina- 
tions of  clamped 
and  simply  supported 

1.4 

Paratlelogrammic 

Simply  supported 

2 

Elliptic 

Clamped 

1.3.4 

Right  triangular 

Clamped  and  simply 
supported 

1.4 

Insofar  as  the  development  of  relatively  new  methods 
of  analysis  is  concerned,  Dharmarajan  and  Chou  (8) 
extended  the  method  of  constant  deflection  lines  - 
originated  for  deflection  [9]  and  free  vibration  of 
isotropic  plates  (10)  - to  orthotropic  plates.  They 
applied  it  to  both  clamped-edge  and  simply-sup- 
portsd-edge  elliptic-pianform  plates.  The  main  dis- 
advantage of  this  method  is  that  it  can  be  used  to 
determine  only  the  fundamental  frequency. 

Another  method,  not  actually  new  but  seldom 
applied  to  vibration  of  plates,  is  the  Kantorovich 
method  (11),  known  as  the  L$vy  method  in  static 
plate  deflection  theory  (12).  This  method  has 
recently  been  applied  to  free  vibration  of  ortho- 
tropic plates  (13) . Solutions  were  obtained  by  finite 
differences  with  respect  to  one  position  coordinate. 
Uniform-thickness  rectangular  plates  with  seven 
different  boundary  conditions  and  linearly  tapered 
rectangular  plates  with  simply  supported  edges 
were  studied.  Values  for  the  first  six  natural  frequen- 
cies were  In  fairly  close  agreement  with  those  ob- 
tained using  other  methods. 

Vljayakumar  and  Remaiah  (14)  used  a variation 
(161  of  Bolotin's  asymptotic  method  (15)  to  obtain 
an  initial  solution;  this  solution  became  the  trial 
function  In  the  Rayleigh  and  Rayleigh-Ritz  methods. 
The  results  for  clamped  rectangular  platas  were  com- 
pared with  those  obtained  previously.  Compared 
with  36-term  Rayleigh-Ritz  results  (17),  the  modi- 
fied Bolotin  estimates  for  the  fundamental  frequency 
were  as  much  as  4.8%  In  error  (always  on  the  low 


side);  the  results  with  the  variation  (14)  are  only 
0.32%  in  error  (sometimes  lower,  other  times  higher). 
It  is  significant  that  a well-known  simple  formula 
(18)  gave  frequency  values  as  much  as  50%  too  high. 

In  addition  to  work  mentioned  above  [13],  varying 
thickness  plates  have  also  been  investigated;  approxi- 
mate fundamental-frequency  expressions  were  de- 
rived for  simply  supported  (19)  and  clamped  (20) 
edge  conditions. 

Sakata  (21)  treated  a multi-bay  continuous  rect- 
angular panel.  He  considered  simple  supports  at  the 
y»0  and  y«b  edges,  free  or  elastically  restrained 
edges  at  x-0  and  x»a,  and  simple  supports  at  the 
Intermediate  supports  x»a/N,  2a/N,  ..  ,(N-1)a/N. 

Laura  and  Luisoni  (22)  used  the  Rayleigh-Ritz 
method  in  conjunction  with  polynomial  modal 
functions  to  analyze  a rectangular  plate  with  dif- 
ferent values  of  elastic  restraint  at  the  edges  and  sub- 
jected to  In-plane  forces.  Beam  functions  were  used 
in  conjunction  with  the  Rayleigh-Ritz  method  to 
treat  in-plane-loaded  rectangular  plates  with  various 
combinations  of  boundary  conditions  (23).  A 
conformal-mapping  technique  used  to  analyze  Ir- 
regularly shaped  plates  under  in-plane  forces  (24) 
led  to  numerical  results  for  the  case  of  a rectangular- 
ly orthotropic  plate  of  circular  planform.  However, 
the  results  appear  to  contain  some  numerical  error 
(25). 

With  the  growing  industrial  importance  of  the  finite- 
element  method,  it  is  not  surprising  that  its  use  in 
composite-plate  vibration  analyses  is  increasing. 
A simplified  mixed  rectangular  element  has  been 
applied  to  free  vibration  of  rectangular  plates  (26) . 
A combination  of  triangular  and  rectangular  elements 
have  been  used  for  the  vibration  of  a cantilever  plate 
representative  of  a delta-planform  airplane  wing  (27) . 

Forced  Vibration 

Relatively  few  investigations  on  forced  vibration  havo 
been  reported  recently,  except  for  some  analytical 
work  (28)  and  experimental  work  (29)  on  respnse 
to  random  excitation. 

An  equivalent  viscous  damping  approach  was  used 
to  predict  response  of  a simply-supported  rectangular 
plate  to  a sinusoidal  point  load  normal  to  the  plate 
at  its  center  (30) . 
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Two  investigation*  concerned  with  thermally  in- 
duced vibrations  were  on  rectangular  (31)  end 
equilateral  triangular  [32]  platei  The  latter  is  severe- 
ly limited  by  the  assumption  that  2(Du  + 2D««)  - 
(Du  Du)'/V  Related  work  on  the  effect  of  thermal 
stresses  induced  by  steady-state  sinusoidal  loading 
has  been  reported  [33] . 

SMALL-DEFLECTION  MOTION  OF  THIN, 
SYMMETRICALLY  LAMINATED. 

RECTANGULARLY  ANISOTROPIC  PANELS 

Relatively  very  little  work  has  been  conducted 
recently.  The  governing  equation  of  motion  for  this 
class  of  panel  is  identical  to  that  in  equation  (1)  with 
the  addition  of  two  bending-twisting  coupling  terms 
on  the  left  side:  4D|  »w,xxxy  and  4D2«w,Xyyy. 

In  optimization  studies  of  rectangular  plates  having 
the  symmetric  balanced  angle-ply  lamination  ar- 
rangement, for  instance  91-91-919,  both  simply 
supported  [34]  and  clamped  [35]  edges  were  con- 
sidered. In  each  case  the  optimal  lamination  angle 
(9)  depends  upon  the  plate  aspect  ratio  as  well  as 
the  laminate  material. 

Laura  and  Grossi  [361  considered  rectangular  plates 
with  supported  edges  having  elastic  rotational  con- 
straints However,  the  latter  conditions  are  satisfied 
only  approximately.  Solution  was  obtained  by  the 
Rayleigh-Ritz  method  with  polynomial  modal  func- 
tions For  cases  in  which  comparisons  with  previous 
solutions  were  possible,  sufficiently  good  agreement 
for  engineering  purposes  was  obtained. 

In  many  aerospace,  marine,  and  automotive  applica- 
tions, cutouts  are  necessary  for  lightening  and  access. 
Thus,  work  on  rectangular  panels  with  rectangular 
cutouts  is  of  practical  importance  [37,  38] . The  cut- 
out is  assumed  equivalent  to  a certain  displacement- 
dependent  external  loading.  Both  simply-supported 
(37)  and  clamped  [38]  edges  and  sinusoidally 
forced  and  free  vibrations  were  treated.  Results 
in  graphical  form  depicted  the  effects  of  plate  ma- 
terial, orientation,  and  plate  and  cutout  geometry. 

In  analyses  of  the  flutter  of  rectangular  panels  under 
high-Mach-number  supersonic  conditions,  numerical 
rasults  were  presented  for  boron-epoxy,  graphite- 
epoxy,  and  boron-aluminum  unidirectional  materials 
at  various  orientations  [39] . 


SMALL-DEFLECTION  MOTION  OF  THIN. 
UNSYMMETRICALLY  LAMINATED. 
RECTANGULARLY  ORTHOTROPIC 
AND  ANISOTROPIC  PANELS 

Finite-element  analyses  are  now  being  used  for  this 
class  of  plate  [40,  41) . The  application  of  the  NAS- 
TRAN  CQDPLT  element  to  rectangular-planform 
cantilever  plates  gave  excellent  agreement  with 
experimental  results  for  boron-epoxy  panels  [40] . 
For  this  set  of  boundary  conditions,  the  reduced 
stiffness  simplification  was  shown  to  be  adequate. 
In  this  approach  the  direct  bending-stretching  coup- 
ling terms  are  omitted,  but  the  actual  bending  stiff- 
ness matrix  [Djj]  is  replaced  by  the  reduced  one 
(Djj*)  defined  as  follows 

[Dij*]  - (Djj]  - [Bjj]  [AijJ-'lBij] 

NASTRAN  CTRIA  2 was  used  to  analyze  the  free 
vibration  of  cantilever  compressor  blades  [41], 
Results  were  in  good  agreement  with  resonant  fre- 
quencies and  nodal  patterns  obtained  experimentally 
by  holography. 

The  Rayleigh-Ritz  method  was  used  to  determine 
the  resonant  frequencies  of  clamped  cross-ply  lamin- 
ated plates  of  circular  planform  [42] . 

The  Galerkin  method  was  used  to  obtain  Fourier- 
form  solution  for  the  flutter  of  arbitrarily  laminated 
thin  plates  [43].  Graphical  results  showed  the  effects 
of  lamination  angle,  lamination  arrangement,  and 
flow  direction  on  the  flutter  parameter  for  plates 
laminated  of  glass-epoxy,  boron-epoxy,  and  graphito- 
epoxy. 


SMALL-DEFLECTION  MOTION  OF  THIN. 
SYMMETRICALLY  OR  UNSYMMETRICALLY 
LAMINATED.  CYLINDR1CALLY  ORTHOTROPIC 
OR  ANISOTROPIC  PANELS 

Apparently  It  would  be  difficult  to  manufacture  a 
uniform-thickness  plate  that  is  cylindrically  ortho- 
tropic; i.e.,  having  the  material-symmetry  axes 
oriented  in  the  radial  and  circumferential  directions. 
Nevertheless,  vibrational  analyses  of  such  plates, 
either  single-iayer/symmetrically  laminated  [44-46] 
or  arbitrarily  laminated  [47]  have  been  reported. 
It  would  seem  to  be  even  more  difficult  to  construct 
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a cylindrically  anisotropic  plat*;  i.e.,  hawing  a full 
array  of  alattic  propartiaa  (no  aymmatriaa  or  zero 
values}  with  respect  to  a cylindrical  coordinate  sys- 
tem. However,  the  vibrations  of  this  class  of  plate 
have  also  bean  analysed  [49] . 

Prat  nap  and  Varadan  (44)  considered  the  axisym- 
metric  free  vibration  of  a solid  circular  plate  of 
cylindrical^  orthotropic  matadal.  Comparison  of  the 
results  obtained  by  the  Rayleigh  energy  method 
(called  the  Lagranglan  method  by  the  authors)  with 
those  obtained  by  the  Galerkin  method  revealed  the 
source  of  some  discrepancies  in  some  existing  analy- 
ses of  this  problem.  The  method  of  Frobenius  for 
the  nonaxisymmetric  case  (45)  was  used  to  obtain 
results  that  have  been  criticized  (44) . 

Ramalah  and  Vijakumar  (46)  extended  previous 
work  to  the  vibration  of  concentric  annular  cylin- 
drical^ orthotropic  plates. 

For  cylindrical^  anisotropic  plates  with  small  flexur- 
al rigidities  and  subject  to  high  in-plane  loads,  the 
numericel-perturbatlon  method  of  matched  asymp- 
totic expansions  was  used  to  study  the  free  vibrations 
(48). 


EFFECTS  OF  THICKNESS-SHEAR 
DEFORMATION  ON  SMALL-DEFLECTION 
MOTION 

It  has  long  been  accepted  that  thickness-shear  defor- 
mation plays  a prominent  role  in  sandwich  panels; 
l.e.,  those  with  one  or  more  thick,  highly  flexible 
cores  and  one  or  more  thin,  relatively  rigid  facings. 
However,  this  effect  Is  often  important  in  composite- 
material  laminated  panels  as  well;  i.e.;  those  with 
all  layers  of  equal  thickness.  In  fact,  it  is  Important 
even  in  single-layer  panels  of  composite  material 
due  to  their  very  low  ratio  of  shear  modulus  to 
Young's  modulus  as  compared  to  Isotropic  materials. 

Symmetric  Laminates 

A Galerk  in-type  approach  was  used  to  study  the  free 
vibration  of  rectangular  plates  of  rectangularly 
ortho  tropic  material  and  subjected  to  any  combina- 
tion of  simply  supported,  elastically  supported,  or 
clamped  edge  conditions  (49).  Rotatory  Inertia 
was  Included. 


Hinton  extended  the  finite-strip  method  to  the 
symmetrically  laminated  anisotropic  case  (50)  and 
also  developed  an  analogous  finite  element  (51). 

The  Mindlin-Goodman  procedure  (53)  was  used  to 
obtain  a solution  for  the  transient  response  of  an 
infinite  long  plate  strip  simply  supported  along  tha 
sides  (52) . The  effects  of  pulse  shape  (rectangular, 
triangular,  and  sinusoidal)  and  pulse  dwell  time  on 
dynamic  load  factor,  maximum  deflection,  maximum 
bending  stress,  and  maximum  Interlaminar  shear 
stress  were  investigated.  A somewhat  similar  analysis 
has  been  reported  for  the  plane-strain  case  in  the 
presence  of  residual  thermal  stresses  (54) . The  results 
indicated  that  thermal  stresses  are  usually  detri- 
mental to  dynamic  behavior. 


Unsymmetric  Laminates 

There  has  been  considerable  activity  in  the  case  of 
free  vibration  of  rectangular-planform  plates.  Bert 
and  Chen  (55)  presented  a closed-form  solution  for 
a certain  kind  of  simply  supported  edges.  The  results 
of  a mixed  finite  element  to  this  problem  (56)  were 
in  good  agreement.  An  eight-node,  40-degreeof- 
freedom  thick-plate  element  was  in  good  agreement 
with  experimentally  determined  natural  frequencies 
(57).  The  British  Program  VIPASA  has  been  ex- 
tended (58) . 

Venkatesan  and  Kunukkasseril  (59)  analyzed  the 
free  vibration  of  a circular -planform  plate  laminated 
of  different  isotropic  materials.  Symbolic  mani- 
pulation techniques  in  conjunction  with  a Rayleigh- 
Rltz  analysis  of  a clamped  elliptic  plate  hare  been 
used  (60,61). 

A thick  multilayer  laminate  theory  has  been  for- 
mulated and  used  to  predict  the  natural  frequencies 
of  a simply-supported  rectangular  plate  (62) . Bend- 
ing, traction,  and  shear  effects  are  Included  for  each 
layer,  and  continuity  of  stresses  and  displacements 
at  the  interlaminar  interfaces  is  maintained.  Dong  and 
Pauley  (63)  used  elements  through  the  thickness 
to  Investigate  plane  wave  propagation  in  thick, 
laminated  plates. 

The  forced  response  of  a symmetric  crose-ply  sub- 
strate was  analyzed  with  a constrained  damping 
layer  (64) . 


Damping  hat  been  incorporated  into  dynamic  analy 
»i  of  sandwich  panels  in  three  ways  In  recent  work: 


Earlier  literature  on  dynamics  of  sandwich  panels 
is  extensive  (1,  05) ; recent  analyses  have  thus  been 
directed  toward  specialised  ejects  - for  example, 
the  effect  of  thteknasanormal  flexibility  of  core 
material  (06).  The  effect  on  impact  is  somewhat 
analogous  to  the  impact  of  a thickness-normal-rigid 
sandwich  plate  on  a free  surface  of  liquid. 


e Um  of  single  composite  Ion  lector  (72,  73) . 
The  main  disadvantage  of  this  approach  is  that 
the  composite  loss  factor  is  most  conveniently 
determined  empirically  from  experiments. 


a Um  of  a complex -modulus  approach  for  the 
core-material  drear  moduli  only  (74,  75) . The 
main  disadvantage  of  this  approach  is  that  it 
categorically  neglects  any  damping  in  the 
facings. 


A rectangular  sandwich  panel  with  unsymmetric 
orthotropic  facings  was  subjected  to  a time-harmonic 
point  loading  normal  to  the  plate  (67) . Calculations 
tor  a sandwich  with  one  facing  of  aluminum,  the 
other  of  steal,  and  a polyvinyl  chloride  core  indicated 
that  the  affect  of  symmetric  modes  on  antisymmetric 
ones  was  not  pronounced. 


e Um  of  a complex-modulus  approach  tor  the 
Young's  modulus  of  the  facing  materiel  as  well 
as  the  sheer  modulus  of  the  core  material  (76) . 
In  any  caM.  this  is  the  most  accurate  model, 
especially  with  facings  of  polymer-matrix  com- 
posite material.  Unfortunately,  however,  the 
work  was  limited  to  isotropic  facings  and  an 
isotropic  core. 


L«e  end  Chang  (68)  investigated  the  dispersion 
relations  for  a sandwich  with  symmetric  facings  and 
isotropic  materials  throughout.  The  material  arrange- 
ment was  typical  of  afoctro  static  ally  charged  pre- 
cipitator plates;  i.e.,  a heavy  .stiff  core  (the  plate 
itMlf)  and  a light,  flexible  facing  (coatings  of  dust 
particles  adhering  to  the  plate).  This  is  the  convene 
of  the  material  arrangement  typical  of  structural 
applications. 


Flutter  of  a sandwich  panel  was  analysed  (77) 
using  the  finite-element  method. 


Relatively  compact  electronic  devices  are  sometimes 
mounted  on  Mndwich  panels.  It  Is  of  interest,  there- 
fore. to  study  the  effect  of  a mass-spring-dashpot 
system  attached  et  an  arbitrary  point.  The  damped 
free  vibration  of  one  such  system  had  a rectangular 
plate  with  symmetric  orthotropic  facings  and  ortho- 
tropic  core  (69) . 


EFFECTS  OF  GEOMETRIC 
NONLINEARITY 


Of  the  two  sources  of  nonlinearity  in  the  dynamics 
of  composite-material  plates,  the  most  common  is 
geometric  nonlinearity  at  finite  deflections  due  to 
the  stiffening  membrane  action  when  the  plate 
edges  are  prevented  from  undergoing  any  in-plane 
diqslecement.  The  effect  is  e hardening-spring  action 
that  cauMS  tht  natural  frequencies  to  incretM  as 
the  amplitude  of  motion  is  increased. 


Chen  and  Came  (70)  conducted  a finite-element  and 
experimental  investigation  of  an  open  Mndwich  panel 
consisting  of  a single  flat  plate  with  a traperoidal 
corrugated  member  attached.  An  epoxy-adhesive 
bonded  structure  had  a greater  stiffness  then  a 
spot-welded  one. 


Thin  Laminates 

A variety  of  solutions  tor  uniform-thickness,  rect- 
angularly orthotropic  thin  plates  undergoing  large- 
amplitude  vibration  have  recently  appeared.  It  is 
gratifying  that  only  one  (78)  of  theM  investigations 
usad  the  overworked  Berger’s  hypothesis,  which  has 
been  shown  to  be  unreliable  in  many  cases  (1). 
Various  methods  have  been  applied  to  rectangular 
orthotropic  plates,  the  Galerkin  technique  (79,  80) . 
the  perturbation  method  (81),  and  e simplified 
1 6-dsgr*e<jf -freedom  rectangular  conforming  finite 
element  (82) . The  letter  kind  of  element  hM  also 
been  usad  with  cylindrical  orthotropic  plates  (83) . 


A number  of  recent  papers  have  been  concerned  with 
the  damping  behavior  of  sandwich  panels.  In  an 
experimental  Investigation  of  a foam -core  panel, 
and  on  the  basis  of  a beam  analysis,  the  system 
wss  believed  to  behave  as  a distributed-parameter 
tuned  dynamic  absorber  (71).  At  low  frequencies, 
the  behavior  was  dominated  by  thickness-normal 
action  analogous  to  a spring-mass  system;  at  higher 
frequencies  the  system  behaved  more  like  a homo- 
geneous beam. 
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Radially  tapered  cylindrieally  orthotropic  circular 
Plata*  war*  treated  analytically  (84)  and  with  finita 
clamant*  (85),  Annular  platat  with  two  different 
nonlinear  taper*  ware  uaad  to  obtain  tolution*  by  a 
combination  of  the  Kantorovich  method  and  the 
Newton- Raphton  iteration  schema  (84).  A linear- 
tapered  solid  circular  plate  was  treated  with  finita 
element*  (85).  Two  large-emplitude  analyses  of 
unsymmetric  laminates  have  appeared  (86,  87). 
In  one  case  antisymmetric  croas-ply  rectangular 
plates  with  two  opposite  edges  simply  supported 
and  the  other  two  edges  clamped  were  studied  using 
the  G dark  in  technique.  Both  stress-free  and  immov- 
able in-plane  edge  conditions  ware  included . 

Most  of  the  analyses  mentioned  above  involved  only 
on*  term  in  position  and  thus  neglected  the  effects 
of  model  coupling.  An  exception  is  the  work  of 
Chia  and  Prabhakara  (87) , who  presented  multiple- 
mode solutions  for  rectangular  plates  of  both  anti- 
symmetric angle-ply  and  antisymmetric  cross-ply 
with  both  simply  supported  and  damped  edges. 
Although  the  effect  of  modal  coupling  on  the  non- 
linear frequencies  of  isotropic  plates  is  not  significant, 
It  can  be  significant  for  composite  panels,  especially 
ctamped-edge  high-modulus  laminates. 

Thick  Laminates  and  Sandwich  Panels 
Sathyemoorthy  (88.  89)  included  the  effects  of 
thickness  sheer  deformation  and  rotatory  inertia  for 
symmetrically  laminated  rectangularly  orthotropic 
plates.  Unfortunately,  he  used  the  Berger  hypo- 
thesis. rather  than  working  directly  with  the  dynamic 
version  of  the  von  Kerman  equations.  The  Berger 
hypothesis  was  also  used  to  derive  a relatively  simple 
frequency  expression  for  a simply  supported  rect- 
angular plate  of  isotropic  symmetric  sandwich  con- 
struction (90) . 

In  an  analysis  of  the  nonlinear  dynamic  sn«>-through 
symmetric  buckling  of  a simply  supported  rectangular 
plate  of  isotropic  symmetric  sandwich  structure,  the 
plat*  Is  assumed  to  have  a slight  initial  curvature 
(91] . A double  Fourier  series  in  position  and  Hou- 
bolt's  timewise  Integration  scheme  was  used. 


EFFECTS  OF  NONLINEAR  MATERIAL 
BEHAVIOR 

Physical  (material  behavior)  nonlinearity  may  be 
due  to  nonlinear  stress-strain  relations  such  as  exem- 


plified by  ordinary  yielding  or  that  due  to  the  soften- 
ing nonlinearity  encountered  in  the  shear  stress  vs. 
shear  strain  relation  for  certain  composites. 

Zak  (92)  considered  physical  nonlinearity  in  the 
form  of  orthotropic  elestoplestic  material  behavior. 
Me  also  considered  lamination  effects,  thickness  shear 
deformation,  and  geometric  nonlinearity.  His  analy- 
sis is  baatd  on  the  use  of  numerical  timewise  integra- 
tion and  a uniform-thickness  quadrilateral  finite 
element  stacked  in  the  thickness  direction  to  repre- 
sent the  composite  layers. 

Sun  and  Shafey  (93)  considered  material-behavior 
nonlinearity  in  the  form  of  an  additional  cubic 
shear-strain  term  in  the  thickness  shear  stress-strain 
relation.  They  investigated  harmonic  wave  propaga- 
tion, formation  of  shock  waves,  and  free  vibrations. 
A much  more  elaborate  analysis  included  up  to 
fourth-order  elastic  constants  - i.e.,  up  to  cubics  of 
the  strains  (94) . 


TRENDS  AND  SUGGESTIONS 
FOR  FUTURE  RESEARCH 

The  following  trends  are  notable  since  the  first  review 
of  this  series  ( 1) : 

* Use  of  the  finite  element  method  for  both 
linear  and  nonlinear  problems  (95) 

e Increased  study  of  such  practical  complications 
as  cutouts,  attached  mass-spring  damper  sys- 
tems, elastically  supported  edges,  and  random 
excitation 

* Increased  activity  in  the  area  of  geometrically 
nonlinear  problems 

The  author  believes  that  the  following  aspects  should 
be  Investigated  more  fully: 

* Mora  general  lamination  schemes  with  attention 
toward  their  use  in  optimal  design  syntheses 
for  such  multiple  conditions  a*  flight  loadings, 
noise,  buckling,  and  flutter 

* More  realistic  material  models,  with  respect 
to  both  effects  of  frequency  and  temperature 
on  stiffnesses  and  damping  and  to  fatigue 


I 

i 


. 


f 

i 


t 

t 


properties  end  even  the  nonlineerlty  of  the 
stress-strain  relations 

e Interactions  between  vibration  and  material 
flaws,  for  example,  crack  propagation  under 
various  kinds  of  vibratory  excitations 

e Experimental  verification  of  analyses,  as  well 
as  vibration  testing  as  a system  identification 
technique 
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